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Abstract 



I discuss a simple numerical algorithm for the direct evaluation of multiple 
Grassmann integrals. The approach is exact, suffers no Fermion sign prob- 
lems, and allows arbitrarily complicated interactions. Memory requirements 
grow exponentially with the interaction range and the transverse size of the 
system. Low dimensional systems of order a thousand Grassmann variables 
can be evaluated on a workstation. The technique is illustrated with a spinless 
fermion hopping along a one dimensional chain. 
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In path integral formulations of quantum field theory, fermions are treated via integrals 
over anti-commuting Grassmann variables Q. This gives an elegant framework for the for- 
mal establishment of Feynman perturbation theory. For non-perturbative approaches, such 
as Monte Carlo studies with a discrete lattice regulator, these variables are more prob- 
lematic. Essentially all such approaches formally integrate the fermionic fields in terms of 
determinants depending only on the bosonic degrees of freedom. Further manipulations give 
rise to the algorithms which dominate current lattice gauge simulations. 

Frequently, however, this approach has serious shortcomings. In particular, when a 
background fermion density is desired, as for baryon rich regions of heavy ion scattering, 
these determinants are not positive, making Monte Carlo evaluations tedious on any but 
the smallest systems f2[|. This problem also appears in studies of many electron systems, 
particularly when doped away from half filling. 

In this note I explore an alternative possibility of directly evaluating the fermionic inte- 
grals, doing the necessary combinatorics on a computer. This is inevitably a rather tedious 
task, with the expected effort growing exponentially with volume. Nevertheless, in the pres- 
ence of the sign problem, all other known algorithms are also exponential. My main result 
is that this growth can be controlled to a transverse section of the system. I illustrate the 
technique with low dimensional systems involving of order a thousand Grassmann variables. 

I begin with a set of n anti-commuting Grassmann variables {ipi}, satisfying = 
ipiipj + ipjipi = 0. Integration is uniquely determined up to an overall normalization by 
requiring linearity and "translation" invariance / dipf{ip + ip') = J dipf(ip). For a single 
variable, I normalize things so that 

[dip V = 1 / dip 1 = 0. (1) 

In summary, to obtain a non-vanishing contribution, every integration variable must appear 
exactly once in the expansion of the integrand. 

Consider an arbitrary action S(ip), a function of these anti-commuting variables, inserted 
into a path integral. In particular, I want to evaluate 

Z = J ' dip n ...dip x e 5W . (2) 

Formally this requires expanding the exponent and keeping those terms containing exactly 
one factor of each ip. 

I first convert the required expansion into operator manipulations in a Fock space. In- 
troduce a fermionic creation-annihilation pair for each fermionic field, i.e. ipi <-> {aj,a«}. 
These satisfy the usual relations 

[oi,a]] + = Sij. (3) 

The space is built up by applying creation operators to the vacuum, which satisfies Oj|0) = 0. 
It is convenient to introduce the completely occupied "full" state 

\F)=al..ata\\0). (4) 

With these definitions, I rewrite my basic path integral as the matrix element 
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Z = (0| e s{a) \F). (5) 

Expanding e s , a non- vanishing contribution requires one factor of eij for each Fermion. This 
is the same rule as for Grassmann integration. 

I now manipulate this expression towards a sequential evaluation. Select a single variable 
ipi and define Si (a) as all terms from the action involving a factor of aj. I define the 
complement Si to be all other terms, so that S = Si + SV For simplicity I assume that my 
action is bosonic so that Si and Si commute. Thus I have 

Z = (0\e St e Sl \F). (6) 

Observe that since Si contains no factors of a^, the occupation number for that variable, 
Hi = ojoi, would vanish if inserted between the two operators in Eq. ([]). I thus can insert 
a projection operator 1 — rii between these factors 

Z=(0\e St (l-O e s >\F). (7) 

The next two steps are not essential to the algorithm, but simplify the appearance of the 
final result. First I use af = and the fact that Si is linear in a^. Thus, the right hand 
factor expands as only two terms 

Z=(0\e Si (1-rO (8) 

Since 1 — rij projects out an empty state at location i, I trivially have (1 — nj) = 0. This 
means Si (1 — rij) = and I can replace Si on the left with the full action 

Z=(0\e s (1-m) (l + Si)\F). (9) 

Finally, I repeat this procedure successively for all variables, giving the main result 







J] ((!-*)(! + Si)) 



(10) 



where I have assumed that the action has any remaining constant pieces removed. 

This summarizes the basic algorithm. Begin by setting up an associative array for storing 
general states of the Fock space. Standard hash table techniques allow rapid storage and 
retrieval. More explicitly, for a given state 

M = £x.l«>. (ii) 

s 

store the numbers \s labeled by the respective states \s). At the outset this table is very 
short, only containing one entry for the full state. The algorithm proceeds with a loop over 
the Grassmann variables. For a given ipi, first apply (1 + Sj) to the stored state. Then empty 
that location with the projector 1 — n{. After all sites are integrated over, only the empty 
state survives, with the desired integral as its coefficient. 

The advantage of the scheme becomes apparent with a local interaction. All sites that 
have previously been visited are empty, and thus involve no information. Any unvisited 
locations outside the range of the interaction are still filled, and thus also involve no storage. 
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All relevant Fock states are nontrivial only for unvisited sites within the interaction range 
of previously visited sites. Sweeping through the system in a direction referred to as "lon- 
gitudinal," we only need keep track of a "transverse" slice of the model. This is illustrated 
in Fig. flU). Thus, although the total number of basis states in the Fock space is two to the 
number of Grassmann variables, the storage requirements only grow as two to the transverse 
volume of the system. 

Note that the algorithm makes no assumptions about the precise form of the interaction. 
The approach is exact, with no sign problems. The complexity does grow severely with 
interaction range, probably limiting practical applications to short range interactions in 
low dimensions. Note that the effort only grows linearly with the longitudinal dimension, 
allowing very long systems in one direction. This discussion has been in the context of "real" 
Grassmann variables. For "complex" variables simply treat ip and ip* independently. 

In the transverse direction the boundary conditions are essentially arbitrary, but longitu- 
dinal boundaries should not be periodic. To make them so requires maintaining information 
on both the top and bottom layers of the growing integration region, squaring the difficulty 
Note that the technique is similar to the finite lattice method used for series expansions || , 
and closely related to attempts to directly enumerate fermionic world lines Q . 

For a simple demonstration, consider a spin-less fermion hopping along a line of sites. I 
introduce a complex Grassmann variable on each site of a two dimensional lattice and study 

(dtfjd^*)e St+Sh (12) 



with the temporal hopping term of form 



and the spatial hopping 



s t = J2^M,t-^,t-i) (is) 

i,t 



S h = kJ2^,A+i,t + ^ +1 , t A,t. (14) 

i,t 



I take N t sites in the time t direction and N spatial sites. Here the one sided form of the 
temporal hopping insures the model has an Hermitian transfer matrix in this direction |J. I 
treat the time direction as my "transverse" coordinate, growing the lattice along the spatial 
chain. After each row with Nt = 10, the number of Fock states stored in the hash table 
rises to 184,756. In Fig. I plot the resulting "free energy" F = \og(Z)/NiN t for k — 1 
on a 50 site chain as a function of the number of time slices. The 50 by 10 case has one 
thousand Grassmann variables. This model is easily solved in the infinite length limit by 
Fourier transform. For infinite Nt this gives 



NN t 



_ 77 

J I dglog(l + 2cos(g)) = 0.388748 . . . (15) 



Note how the results in the figure oscillate about this line, indicating that the transfer 
matrix, while Hermitian, is not positive definite. The transfer matrix with this action 
becomes positive definite only for \k\ < \. 
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Now I make the model somewhat less trivial with a four Fermion interaction. I take 
S = S t + Sh + Sj with 

Si = gJ2 ^tt^i^t+i^i+u- (16) 

i.t 

In a transfer matrix formalism this represents an interaction Hamiltonian of form Hj = 
Yli n i n i+i- Bosonization relates this Hamiltonian to the anisotropic quantum Heisenberg 
model, but this equivalence is not being used here. The above algorithm requires essentially 
the same computer resources as the free case. In Fig. (|^) I show the N t dependence for 
the free energy with k = 1 and g = ±1. To reach N t = 13 for this figure, I reduced 
memory requirements further by using temporal translation invariance after each layer was 
integrated. 

With Monte Carlo studies of many fermion systems, the introduction of a chemical 
potential term can be highly problematic due to cancellations. Here, however, it is just 
another local interaction of negligible cost. As an illustration, to the above model I add 
another term and take S = S t + Sh + Si + Sm with 

S M = M^ttAt- (17) 

i,t 

I can use this term to regulate the "filling," which can be approximately monitored as 
(1 + M)-j£. Here I include the extra factor of 1 + M to compensate partially for finite N t 
artifacts. This quantity is shown as a function of MN t in Fig. (^), obtained on an N t = 8 
by Ni = 20 lattice with a spatial hopping parameter of k = 0.1. For this figure I made a 
crude extrapolation in chain length by defining F — log (-S^zu) ■ Note how the four 
fermion coupling enhances the filling. The crossing of the curves at large chemical potential 
is a consequence of strong coupling at finite N t . 

An obvious system for future study is the Hubbard model 0. This requires 4 Grass- 
mann variables per site corresponding to ip* and ip for spins up and down. Higher spatial 
dimensions strongly increase the size of the transverse volume and will limit practical system 
volumes, but this may be compensated for by the lack of sign problems. 

Another potential application relates to the use of Shockley surface states to formulate 
chiral gauge theories in lattice gauge theory. These so called "domain wall Fermions" have 
unresolved questions due to a natural pairing of surfaces. One suggestion |7j proposes a 
four fermion coupling on one surface to remove the spurious modes. The effects of such a 
coupling might be studied in a truncated model via the above techniques. 
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FIG. 1. Integrating out sequentially, all finished sites are empty and out of range sites are filled. 
When integrating the site labeled ip, only those sites labeled "X" are undetermined. 
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FIG. 2. The free energy F = \og{Z)/NiNt on an Ni = 50 site chain as a function of the number 
of time slices Nt- The infinite volume solution is shown by the horizontal line. 
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FIG. 3. The free energy F = log(Z) /NiNt with a four fermion interaction as described in the 
text, plotted as a function of the number of time slices Nf. The chain has iVj = 50 sites. Points 
are shown for k = 1 and g = ±1. 
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FIG. 4. The occupancy of a 20 site chain as a function of the chemical potential scaled by the 
number of time slices. Note how the filling occurs earlier or later depending on the sign of the four 
fermion coupling. Here the spatial hopping parameter is taken as K = 0.1 
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